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Abstract: A family of cosmological solutions with Ricci-flat spaces in the theory with 

several scalar fields and multiple exponential potential is obtained when coupling vectors in 
exponents obey certain relations. Two subclasses of solutions with power-law and exponential 
behaviour of scale factors are singled out. It is proved that power-law solutions may take place 
only when coupling vectors are linearly independent and exponential dependence occurs for 
linearly dependent set of coupling vectors. A subfamily of solutions with accelerated expansion 
is singled out. A generalized isotropization behaviours of certain classes of general solutions 
are found. In quantum case exact solutions to Wheeler-DeWitt equation are obtained and 
special "ground state" wave functions are considered. 
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1. Introduction 

Recently, the discovery of the cosmic acceleration [1, 2] stimulated a lot of papers on multidi- 
mensional cosmology with the aim to explain this phenomenon using certain multidimensional 
models, e.g. those of superstring or supergravity origin (see [3]- [12], and references therein). 
These solutions usually deal with time-dependent scale factors of internal spaces, thus over- 
coming a "no-go" theorem for static (and compact) internal spaces [13]. It should be noted 
that certain part of publications does not deal with new exact solutions but use old ones 
(sometimes rediscovered or written in different parametrization). 
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A simple example of vacuum solution with acceleration was considered by Townsend and 
Wohlfarth in [3] . This is a (4 + n)-dimensional solution to vacuum Einstein equations with 
n-dimensional internal space of negative curvature and 4-dimensional "our space" containing 
expanding 3-dimensional flat subspace. The solution from [3] is a special solution, in a special 
model with a very special choice of conformal frame. We note that it is a special case of more 
general vacuum solution from [14], describing the "evolution" of (n — 1) Ricci-flat spaces and 
one Einstein space of non-zero curvature. The solution from [14] was generalized to scalar- 
vacuum case in [15, 16, 17] and to the case of composite p-brane configurations in [19] (for 
non-composite case see also [18]). For a review see also [20]. We note that the solutions 
from [14, 16] may be also used for generation of special solutions with several curved factor 
spaces (e.g. from [43]) using a curvature-splitting trick [35] (when Einstein space of non-zero 
curvature is chosen as a product of several Einstein spaces). (For three integrable classes of 
vacuum cosmological solutions with two factor spaces of non-zero curvature in dimensions 
D = 10,11 see [35].) 

At present rather popular models are those with multiple exponential potential of the 
scalar fields (see, for example, [21, 23, 24, 25, 41, 42, 44] and refs. therein). 

Such potentials also arise naturally in certain supergravitational models [22] and in sigma- 
models [26], related to configurations with p-branes. 

Here we consider the D-dimensional model governed by the action 

S act = [ d D z v ^\{K[g]-h a pg MN dM<p a d N <pP-2V ip (<p)} + S GH , (1.1) 
with the scalar potential 

V M = ^A s exp[2A s (v?)]. 

Here 5gh is the standard Gibbons-Hawking boundary term [36] 
The notations used here are the following ones: 

* if = (ip a ) is the vector from scalar fields in the space El' with a metric determined by 
non-degenerate I x I matrix {h a p) with inverse one (h alS ), 

a = l,...,l; 

* A s are constant terms, 
s e S; 

* X s is an 1-form on JR l : \ s (ip) = X sa f a ; A° = /i a/3 A S/ 3; 

* 9 = 9MNdz M dz N is the metric, \g\ = \ det(gMN)\, 
M, N are world indicies 

(that may be numerated by 1,...,D); 
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(1.2) 



* i,j = 0,...,n are indicies describing a chain of factor spaces; 
A = i, a and B = j,(3 are minisuperspace indicies 
(that may be numerated also by 0, . . . , n + I). 

In this article we obtain new (and general) families of classical and quantum "cosmolog- 
ical" solutions with vector coupling constants obeying 

XsaK'ph^ = (1.3) 

for s / s' and 

XsaXs^ + (1.4) 

Here we find exact solutions with scale factors and scalar fields depending upon one 
variable u ("time"). We keep the parameter w = ±1 in the metric g = we 2l< ^du <8> du + . . . 
for a future consideration of static configurations [40]. They may be of interest in a context 
of "black hole" solutions with unusual asymptotics (see [27, 28, 29, 30, 12] and references 
therein). We single out two new subclasses of solutions with power-law and exponential 
behaviour of scale factors, e.g. those with accelerated expansion. A new result here is 
generalized isotropization behaviours of certain classes of classical solutions that give new 
examples of asymptotically accelerated expansion. 

The paper is organized as follows. In Section 2 the model is formulated. In Section 3 a 
class of classical cosmological solutions corresponding to orthogonal [/-vectors is presented. 
Section 4 deals with special solutions exhibiting power-law and exponential dependence of 
scale factors. Here we describe the solutions with accelerated expansion. A special subsection 
is devoted to a generalized isotropization. In section 5 quantum analogues of the classical 
solutions are considered. In Appendix the equations of motion, corresponding to (1.1) are 
written (Appendix A) and derivations of special solutions of power law and exponentional 
type are presented (Appendix B and C); also a classification of vectors in Euclidean space 
obeying relations (1.3) and (1.4) is done (Appendix D). 

We note that the classical and quantum solutions for cosmological constant case (i.e. for 
one term in potential with A s = 0) were considered earlier in [31] and [32] for vacuum and 
scalar-vacuum cases, respectively. 

2. The model 

Let 

M = R. x M x ... x M n (2.1) 
be a manifold equipped with the metric 

n 

g = we 2 ^du0du + Y,^ iu) 9\ ( 2 -2) 
i=0 
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where w = ±1, 1R. is open interval in IR, u is a distinguished coordinate; g % is a Ricci-flat 
metric on cij-dimensional manifold Mf 



di = dimMj, t = 0,...,n;neN. 
For dilatonic scalar fields we put 



K mi n i [9 i ]=0, (2-3) 



<p a = <p a {u), (2.4) 



It may be verified that the equations of motion corresponding to (1.1) for the field con- 
figuration (2.2), (2.4) are equivalent to equations of motion for 1-dimensional a- model with 
the action 



(2.5) 



Sa = \f duX^Ci.jO'o' + K^ a ^ - 2AT 2 V j, 
where x = dx/du, 

V = -wV^e 2 ^ (2.6) 
is the potential (V v is defined in (1.2)) with 

n 

7O (0) = ^>^, (2.7) 

i=0 

and 

AA = exp( 7o (0) -7) >0 (2.8) 

is lapse function. Here 

Gij = - dd 3 , G ij = ^ + , (2.9) 

i, j = 0, . . . ,n, are components of a gravitational part of minisupermetric and its dual [33]. 
2.1 Minisuperspace notations 

In what follows we consider minisuperspace IR n+1+ ' with points 

x = (x A ) = (4>\<p a ) (2.10) 

equipped by minisuperspace metric G = Gab(1x a ® dx B defined by the matrix and inverse 
one as follows: 
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The potential (2.6) reads 

V = -w^A s e 2US ^ (2.12) 
where II s (x) = U\x A are defined as follows 

n 

U s = U S (<P, ip) = \ sa <p a + di<P\ (2-13) 

or in components 

(U s A ) = (d i ,X sa ). (2.14) 

The integrability of the Lagrange system (2.5) depends upon the scalar products of co- 
vectors II s corresponding to G: 

(U,U') = G AB U A U' B , (2.15) 
These products have the following form 

(U s ,U s ') = -b + X s -X s ,, (2.16) 

where s,s' G S and 

A s • X s i = X sa X s > /3h al3 \ (2-17) 

, D-l i \ 

b =D^2- (2 - 18) 

We put the orthogonality restriction on the vectors II s 

(II s , II s ') = -b + X s ■ X s > = 0, (2.19) 
s s', and also consider non-degenerate case 

(U s ,U s ) = -b + X s -X s ^0 (2.20) 

(see (1.3) and (1.4)). 

In what follows we denote 

fc. = (C^ C^" 1 = (2-21) 
The further consideration is based upon the orthogonality conditions assumed. 
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3. Classical solutions 



Here we will integrate the Lagrange equations corresponding to the action (2.5) in the har- 
monic time gauge 7 = 7o- We get a Lagrangian 



1 



L = -G AB x A x B - V, 



(3.1) 



where {Gab) and V are defined in (2.11) and (2.12), respectively. The zero-energy constraint 
reads: 



1 



E = -G AB x A x B + V = 0. 



(3.2) 



The solutions to Lagrange equations corresponding to (3.1) in the orthogonal and non- 
degenerate case (2.19), (2.20) read (see [34, 19]): 



x 



— E 



U s 



s&S 



(U S ,U S 



In |/ s | + c A u + c A , 



Here 



f s = R s sh(VC s (u - u s )), C s > 0, e s < 
R s sm(^/\c7\(u - u a )), C s < 0, e s < 
R s ch(v^(u-u s )), C s >0, e s >0 
\2As/hs\ 1/2 (u-u s ), C s = 0, e s <0 



where u s and C s are constants and 



(3.3) 



(3.4) 
(3.5) 
(3.6) 
(3.7) 



R s = \2A s /(h s Cs)\ 1/2 , 
e s = -sign(wA s h s ). 

Vectors c = (c A ) and c = (c A ) satisfy the linear constraint relations 



(3.8) 
(3.9) 



U s (c) = U s A c A = d iC l + X sa c a = 0, 

i=0 
n 

U s (c) = U%c A = di& + X sa c a = 0. 



i=0 



The zero-energy constraint reads 

E = Y,E S + \g AB c A c b = ^ 



(3.10) 
(3.11) 



(3.12) 
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where C s = 2E S (U S , U s ), or, equivalently, 

n / n \ ^ 

J] + + ^(c 1 ) 2 - £ d t c l = 0. (3.13) 

seS i=0 \i=0 / 

The solutions. 

Using (3.3) and relations for contravariant components U sA = G AB \J S B 

v si = G ij u s = __i_ ; U sa = X^, (3.14) 
we are led to relations for the metric 

SeS ' seS »-n J 



and scalar fields 



^ = -^M"ln|/ s |+c Q « + c a , (3.16) 

a = 1, . . . , I. 
Here 

n n 

c = Y^dic\ c = ^2dic\ (3.17) 

i=0 i=0 

The functions f s are defined in (3.4)-(3.7) and the relations on the parameters of solutions 
c A , c A (A = i,a) and C s imposed in (3.10)-(3.11) and (3.13), respectively. (h s are defined in 
(2.21).) 



4. Special solutions 

Now we consider a special case of classical solutions from the previous section when C s = 
u s = c 1 = c a = and 

wA s (X 2 s -b) >0, (4.1) 

s G S. 

We get two families of solutions written in synchronous-type variable with: 

A) power-law dependence of scale factors for B ^ — 1, 

B) exponential dependence of scale factors for B = — 1, 
where 



B = B(X) = b^h s . 

seS 



(4.2) 



Note that h s = h s (X s ) = (A 2 - 6) _1 . 
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4.1 Power-law solutions 



Let us consider the solution corresponding to the case B / -1. The solution reads (see 
Appendix B) 

n 

g = wdT®dT + T 2vS ^A ig \ (4.3) 

i=0 

1 „ 

a _ 



s&S 



L^MZlnT + yfi, (4.4) 



where r > 0, 



" = (5 + lj(2>-l)' (45) 



and 



|2A a /fc a | exp(2A aa¥ >g) = (5 + 1)"', (4.6) 

s 6 5; and .Aj > are arbitrary constants. 

Solutions with "acceleration". Let us consider the cosmological case w = — 1. We 
get an accelerated expansion of factor spaces if and only if v > 1 or, equivalently, 

- b < B(X) < -1. (4.7) 

Remark 1. For — 1 < B = B(X) < we get v < 0, that corresponds to a contraction 
of factor spaces. When B < —b or B > we are led to the expansion with deceleration, since 
< v < 1 in this case. The case B = corresponds to static toy "universe". For B = — b we 
obtain the expansion with zero acceleration (u = 1). 

Let us consider the case when the matrix {h a p) is positive definite (e.g. h a p = 5 a p). It 
was shown in Appendix D that the condition B(X) / — 1 in this case implies that vectors 
X s are linearly independent. It follows also from Appendix D that the condition B(X) < —I 
implies that 

A^ () < b, X 2 s >b for all s + s . (4.8) 
Here sq G S. These relations combined with inequalities (4.1) lead to 

A So > 0, A s < for all s / s , (4.9) 

i.e. one exponent term in potential should be positive and others should be negative. 
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4.2 Solutions with exponential scale factors 

Here we consider the solution corresponding to the case B = — 1. The solution reads (see 
Appendix C) 

n 

g = wdr ®dr + exp(2mr) ^ A^g\ (4.10) 

i=0 

( p a = (D-l)mrJ2hsK + ^, (4.11) 

where 

\2A s /h s \ exp(2A SQ ^) = m 2 (D - l) 2 . (4.12) 

s £ <S; m is parameter and .Aj > are arbitrary constants. 

In the cosmological case w = -1 we get an accelerated expansion of factor spaces if and 
only if m > 0. 

Let us consider the case when h a p is positive definite. As it was shown in Appendix D 
the condition B(X) = —1 in this case implies that vectors X s are linearly dependent, the linear 
term in (4.11) vanishes and hence 



<P a = <PS, (4-13) 

i.e. for (h a p) of Euclidean signature all scalar fields are constants when exponential 
expansion is under consideration. We note that any subset of vectors X s ,s 6 S \ {so} is 
linearly independent one in this case (see Appendix D). 

Remark 2. It may be shown that for the solution (4.10)-(4.12) the kinetic term for 
scalar fields vanishes, i.e. h a pg MN Om^On^ 13 = and the scalar potential is constant: 
V v (ip) = V v (ipo) = A and (-u>)2A = m 2 (D - 1)(D - 2) (see (4.1)). For w = -1 and A > 
the metric (4.10) is coinciding with that for multidimensional model from [31] with pure 
cosmological term A. Thus, the exponential expansion of factor spaces is driven by effective 
cosmological term A and is not sensitive to time dependence of scalar fields, that takes place 
only for non-Euclidean signatures of (h a p). 

4.3 One-exponent example 

The simplest case with one term in the potential (1.2) is interesting due to the absence of the 
restriction (2.19). It follows from eq. (4.1) that the potential is positive: 

^ = Ae 2A ^", A>0 (4.14) 

(we discard the indices in this example). 

The accelerated power law expansion takes place, if 

< A 2 < -pi^, (4.15) 
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i.e. A 2 should be small enough but not equal to zero. The inequality is coinciding (up to 
notations) with the eq. (1.9) from [22]. In this case (4.5) reads as v = (X 2 (D — 2)) _1 . 
For the exponential expansion eqs. (4.1) and B(X) = —1 imply 

A 2 = 0, (4.16) 

i.e. A should be zero vector for positive definite matrix {h a p) or light-like vector for (h a p) of 
pseudo-Euclidean signature. 

4.4 Generalized isotropization 

In the case of one exponent the special solutions considered above are attractors in the limit 
r — ► +00 (r is "synchronous" variable) for general solutions with sinh-dependent function f s 
(see (3.4)), when e s < 0, or, equivalently, (4.1) is satisfied and B < —1. This fact may be 
readily verified using the relation f s ~ (u — u s ) for u — ► u s . 

Here we consider a more general case of this attractor behavior for multi-exponent case. 
Let S- be a subset of all S satisfying e s < 0, or, equivalently, relations (4.1). Let us study 
a family of solutions with positive C s and hence with sinh-dependence of f s for s £ S_ and 
cosh-dependence of f s for s £ S \ S-. Consider the solutions on the interval (u#, +oo), where 

u* = max(n s ,s £ S-). (4-17) 

We denote by 5* a subset of all s £ S_ obeying u s = u*. All functions f s (u — u s ) are 

smooth on this interval and near they behave as f s ~ (u — ) for s £ S 1 * , or as (non-zero) 
constants for s £ S \ S*. 

Introducing the synchronous time variable by relation 

T = / duex.p(cu + c)Y]_(fs{u-u s )) bh3 + const (4.18) 

^ s 

we get r — > +oo as u — > ti* in two cases: 

A) B* < -1 and 

r~(u-u») B * +1 ; (4.19) 

B) B* = -1 and 

t ~ — In | ii — Ujf | . (4.20) 

Here we denote 

5, = B(\,S*) = h s b. (4.21) 

Case A. Using arguments analogous to those presented in Appendix B we get the fol- 
lowing asymptotical relations in the limit r — > +oo for B* < — 1 

n 

g as = wdr®dT + Y^ A i T2u *9\ (4.22) 

i=0 
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<s = ~b~T~\ E M>r + ^, (4.23) 

where 

Vtf = (4.24) 

and > are constants. 

Thus, we get an isotropization behaviour of scale factors for r — > +oo, that is supported 

by exponential terms in the potential labelled by s 6 5,. Other exponential terms do not 
contribute to the power-law index . 

For 

- b < B* < -1 (4.25) 

we obtain the accelerated expansion for large enough values of r. 
Case B. 

An analogous consideration in the case B* = — 1 leads to the following asymptotical 
relations in the limit r — > +oo (see also Appendix C) 

n 

g as = wdr fStdT + ^Ai exp(2mr)^, (4.26) 

i=0 

^ = (D-l)mr^W: + * (4-27) 



where m > and 



\2A s /h s \ exp(2A sa <^) = m 2 (D - l) 2 . (4.28) 



s£S». Thus, we obtain the accelerated expansion for large enough values of r. 

This is another regime of isotropization behaviour (for one-exponent case with X s = 
see also [31, 32]). 

Remark 3. It may be shown that for the asymptotical solution (4.26)-(4.28) the kinetic 
term for scalar fields asymptotically vanishes, i.e. h a pg MN ' Bm^^n^v ~~ * as T ~* +°° an d 
the scalar potential is asymptotically constant: 

V v {<p) E As ex P[2A s (^o)] = A* (4.29) 

ses, 

as r — > +oo. Here (— u^A* = m 2 (D — 1)(D — 2). For w = —1 ans A* > the "asymp- 
totical" metric (4.10) is coinciding with that for multidimensional model from [31] with pure 
cosmological term A*. When (h a p) has Euclidean signature the exponential isotropization 
behaviour takes place only for S = S* (it follows from Appendix D). In this case ip as = ipo. 

4.5 Kasner-type behaviour 

Let us study the solution with sinh-dependence in the limit u — > oo. We restrict ourselves by 
Kasner-like behaviour of solutions for small r. The functions (3.4) and (4.18) behave as 
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f s ^ e VcSu^ r _ e (;,s+ c )n ) U ^ +0C) (430) 
where S = ^ h s ^fC~ s . To get the limit r — > +0 we put 

6S + c < 0. (4.31) 
The metric and scalar fields in the limit r — > +0 read as follows: 



5as = wdr ® dr + , (4.32) 

i=0 

¥ 4 = a^lnr + ^, (4.33) 
where the Kasner parameters are given by the formulas 

p = S£S (4 } 

One can see that these parameters obey the generalized Kasner relations: 

n n 

dia 1 = 1, + Kp^a 13 = 1. (4.36) 

i=0 ?=o 

We note that for + c < we get a Kasner- like asymptotics in the limit r — > — oo (with 
r replaced by |r| in (4.32) and (4.33)). 

4.5.1 Example: two spaces with one scalar field 

As an example let us consider the model with two factor spaces, one scalar field and one term 
in potential. It is not difficult to verify, that the parameters 

d (d + di)' di(d + di) 

where R = \J d~od~i(D — 2), satisfy the eqs. (4.36). 

For di = 3, d,2 = 6 the behaviour of scale factors are presented in Fig. 1, when: a) 
A 2 = 1/16 (asymptotically power law expansion); b) A = (asymptotically exponential 
expansion) . 

Remark 4. It will be shown in a separate publication that Kasner-like behaviour near 
the singularity (as r — > +0) is a generic one for n > and a wide class of scalar potentials 
(1.2). 
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*i.*~T /// 
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Figure 1: The behaviour of scale factors of two factor spaces of dimensions do — 3 and d\ = 6 in the 
presence of one scalar field and one term in potential for the cases: a) A 2 = 1/16; b) A = 0. In both 
cases the Kasner parameter for the scalar field is chosen to be zero. 

5. Quantum solutions 

Here we find the quantum analogues of the classical solutions presented above. 
5.1 General solutions 

The standard quantization of the energy constraint [iE = 0, where \i ^ is a parameter, leads 
to the Wheeler-DeWitt (WDW) equation (see, for example, [33, 37, 38]) in the harmonic-time 
gauge 7 = 70 (0) 



-— A[G\+fj.V)9 = 0. 



(5.1) 



Here V is potential defined in (2.12), A[G] is the Laplace operator, corresponding to G, 
and \& is wave function. 

The minisuperspace metric may be diagonalized by the linear transformation 



S A B x B , (z A ) = (z»,z a ,z s ) 



as follows 



G A Bdx A (g> dx B = -dz° ® dz° + V] i]sdz s (g) dz s + dz a ® dz b r] ab , 

ses 



(5.2) 
(5.3) 



where a, b = 1, . . . , n + I - \S\; r] ab = 7/ aa <5 ab ; r] aa = ±1, r? s = sign(/i s ), and 



with 



sZ s = U s (x), 



|/ is |- 1 / 2 = ^Af^>0, 



(5.4) 
(5.5) 
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s G S. 

We are seeking the solution to WDW equation (5.1) by the method of the separation of 
variables, i.e. we put 

* = ( II ) (5 ' 6) 

The wave function (5.6) satisfies to WDW equation (5.1) if 

2H S * S = l^-Vs-^ - 2wA s e 2 ^ * s = 2£ s * s , (5.7) 

where 

V ab p aPb + 2Y,£s = 0- (5.8) 

(compare with classical relation (3.12)). 

The linearly independent solutions to eq. (5.7) reads 

* s (z s ) = B^ U-2w7j s A s ^—] , (5.9) 

where 

lj s = ^/-2S s h s , (5.10) 

s £ S and B^ = a s I w + b s K w are superpositions of modified Bessel function. 

The general solution to the WDW equation (5.1) is a superposition of the "separated" 
solutions (5.6). 

5.2 Special solutions 

Let us consider the special "ground state" solutions with p a = £ s = in (5.9), when relations 
(4.1) are imposed. These solutions read 

* = [J (asMVW^s~\ ve XsM ) + b s H( 1 \^h s T s ~\ ve x °^ , (5.11) 

(11 n 
where Jo and Hq ' are the Bessel and Hankel functions, respectively, and v = exp(^] di(j) 1 ) 

i=0 

is a volume scale factor. These quantum solutions correspond to special classical solutions 
considered in the previous section. 
We note that for large values of 



v s = y/2\h s A s \ v exp(A s (v?)) (5.12) 

we get a quasi-classical regime, that may be obtained (along a line as it was done in [39]) 
using asymptotical relations for Bessel functions. For small values of " quasivolumes" v s we 
are led to a quantum regime. The crucial point here is that the quantum domain is not 
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defined only by volume scale factor v, but by " quasivolumes" v s , depending also upon scalar 
fields. In strong enough scalar fields with certain signs (and the direction of vector ip) and/or 
small |/i s A s | one can obtain a quantum behaviour of a wave function for big enough values of 
volume scale factor v. Analogous "effect" takes place also for quantum solutions with "perfect 
fluid" [16] when certain equations of state are adopted. 

6. Conclusions 

Here we obtained a family of multidimensional cosmological solutions with (n + 1) Ricci-flat 
spaces in the theory with several scalar fields and multiple exponential potential. 

The classical and quantum solutions are obtained if the (orthogonality and non-degeneracy) 
relations (2.19) and (2.20) on [/-vectors are imposed, or, equivalently, when relations (1.3) 
and (1.4) on coupling vectors are satisfied. These solutions in fact correspond to A\ + ... + A\ 
Toda-like solutions (for perfect fluid case see [34]). 

Here we singled out the solutions with power-law and exponential behaviours of scale 
factors. We proved that power-law dependence may take place only when coupling vectors 
\ s are linearly independent; exponential dependence may occur only for a linearly dependent 
set of A s obeying the condition B(X) = —1. We obtained the restriction on coupling vec- 
tors (4.7) that cuts power-law solutions with acceleration. In subsection 4.3 the generalized 
isotropization behavior of certain class of general solutions was found. Any asymptotics with 
power-law or exponential behavior corresponds to a certain subset of exponential terms in 
the potential (all other exponential terms do not contribute to this asymptotical behaviour). 

In the quantum case we solved the Wheeler-DeWitt equation and singled out special 
"ground state" solutions that are quantum analogous of special classical solutions from Section 
4. These solutions are defined as products of Bessel functions and depend upon a set of 
quasivolumes that may be small enough (i.e. belonging to "quantum domain") even when 
the volume scale factor of the toy "universe" is in the classical region (in the Planck scale). 
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. Appendix 

A Equations of motion 

Here we outline for the sake of completeness the equations of motions corresponding to the 
action (1.1) 

TImn - t;9mnTI = Tmn, (A-l) 



- 15 - 



A bb a - Yl 2A^e 2As ^A s = 0. (A.2) 

In (A.l) 

T M n = h a /3 (dM<p a dN<pP ~ ^9MNdp(f a d P (p^ - V^qmn- (A.3) 
B Power- law expansion 

Here we derive the solution (4.3)-(4.6) from the general one. We start from the relation (3.3) 
that can be written in our case as follows 

x A (u) = - h s U sA H^\2K s /h s \u) + c A , (B.4) 

s&S 

where u > 0. 

Introducing a new variable r > by formula 

u = Ct 1 / { - b+1 \ C>0, (B.5) 

we rewrite (B.4) in the following manner 



x A 



= -^2h s U sA \nT + x A (B.6) 

ses 



with constants 

x o = dA ~ E h ' U ' A ln (#WMC) . (B.7) 

Due to orthogonality of ?7 s -vectors the constraints on integration constant U s A c A = may be 
written in the equivalent form 



U%x A = - ln(V|2A s /^|C). (B.8) 
In components the solution (B.6) for (x A ) = (</) l ,ip a ) reads 

b lnr + 4, (B.9) 



(B + 1)(D_1) 
1 

"5 + 



<p a = --J— £ M? lnr + y>g. (B.10) 



Here (x A ) = (^,^). 



For 70 (0) = d i4> % we get 

i=0 



7o(0)= CBTT) lnr + 00 ' (R11) 
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n 

where <p Q = di</> l . 

i=0 

From the definition of " synchronous" variable 

exp(2j o (0))du 2 = dr 2 (B.12) 

we obtain 

C = |5 + l|exp(-0 o ) (B.13) 

and hence (B.8) reads 

X sa ^ = - ln(^/\2A s /h s \C). (B.14) 

The solution (4.3)-(4.6) follows from the formulas (B.9), (B.10), (B.12), (B.14) and A* = 
exp(2$,). 

C Exponential expansion 

The solution (4.10)-(4.12) may be obtained just along a line as it was done for the power-law 
case. The only difference here is the relation 

u = Cexp(Mr), C>0 (C.15) 

instead of (B.5). Using a procedure analogous to considered hereabove we get 

+ (C.16) 

<p a = -Y,hsKM T + ^. (C.17) 

ses 



and 



From (B.12) we get 



7o($ = -Mr + 00- (C.18) 



C=|M|- 1 exp(-^ ) (C.19) 
and introducing new parameter m = —M/(D — 1) we obtain the solution (4.10)-(4.12). 

D b— proper set of vectors 

The statements of this subsection are based on relations (1.3) and (1.4)), that are used in 
Definition of 6-proper set (see below). The conditions for exponential expansion follow from 
Theorem. Lemma 2 is useful for analysis of "power-law" expansion. 

In what follows {Aj} means a set of vectors Ai, . . . , A m G 1R/, I € IN. The vector space IR' 
is equipped with a scalar product (it corresponds to the positive definite matrix (h a p) in the 
bulk of the article) . 

Definition. Let b > 0. The set of vectors {Aj} is called a b— proper one if 
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A; • Xj = 6, (D.20) 
A? + b, (D.21) 

Lemma 1. Let K = {\- Xj) be a matrix of scalar products for a 6— proper set of vectors 
{Ai}- Then 

m 

detK={B + l)]J(X 2 i -b), (D.22) 



i=i 



where 



m , 

B = B(X 1 ,...,X m ) = J2jrri- (D-23) 

i=i * 

Proof. Using relations 

i^- = A, • A, = (A? - + 6 = (A? - 6) + ^—^j (D.24) 
we represent the matrix K as a product of a matrix K and a diagonal matrix D: 

^ = S *i + ' Dij = 6ij (A * " b) ■ (D - 25) 

A product of determinants detD and det^T yields (D.22). Here we used the relation 
detK = B + 1, that can be readily proved. 

Theorem. Let {Aj} be a 6— proper set. The set {Aj} is linearly dependent if, and only if 

B = B{X 1 ,...,X m ) = -l. (D.26) 

In this case 

m . 

A,; 



£aF.=° <D - 27) 

i=l 1 

and any subset of m — 1 vectors from {Aj} is linearly independent one. 

Proof. The linear dependence of vectors {Aj} is equivalent to deti^T = that is equivalent 
(due to Lemma 1) to relation (D.26). 
From the relation 

B(Ai, A m ) = B(Ai, A m _i) + = -1 (D.28) 

A m - 

we get 5(Ai, A m _i) ^ —1 and hence the vectors Ai,...,A m _i are linearly independent. It 
is obvious that any m — 1 A-vectors are also linearly independent. 
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To prove (D.27) let us denote 



rn 



X 2 -b' 



(D.29) 



i=i 



Using the scalar products (D.24) we obtain 



A 2 



B(B + l)/b = 



(D.30) 



and, hence, A = 0. The Theorem is proved. 

As consequence we obtain the inequality on the number of 6-proper vectors: m < I + 1. 
(The equality m = I + 1 takes place when the vectors are linearly dependent.) 

Lemma 2. In the b— proper set {Aj}: 

a. there are no two or more vectors Aj obeying A? < b; 

b. B > — 1 44> vectors Aj are linearly independent and all Xf > b. 

Proof, a. Let us suppose that there are two vectors Aj, Xj, such that A 2 < b and A 2 < b. Then 
b 2 = (Aj • Aj) 2 < A 2 A 2 < b 2 , i.e. we come to a contradiction. 

b. If B > — 1 then by the Theorem the vectors are linearly independent and hence the 
matrix (Aj • Xj) = (Kij) is positive definite. This implies detK > 0. Due to relation (D.22) 
and part a. of this lemma we get: A? > b for all vectors. 

Now, let the vectors {Aj} be linearly independent and all X 2 > b. Then det K > 0, and 
due to (D.22), B + 1 > 0. The lemma is proved. 
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